Vehicular traffic flow at an intersection with the possibility of turning 
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We have developed a Nagel-Schreckenberg cellular automata model for describing of vehicular 
traffic flow at a single intersection. A set of traffic lights operating in fixed-time scheme controls 
the traffic flow. Open boundary condition is applied to the streets each of which conduct a uni- 
directional flow. Streets are single-lane and cars can turn upon reaching to the intersection with 
prescribed probabilities. Extensive Monte Carlo simulations are carried out to flnd the model flow 
characteristics. In particular, we investigate the flows dependence on the signalisation parameters, 
turning probabilities and input rates. It is shown that for each set of parameters, there exist a 
plateau region inside which the total outflow from the intersection remains almost constant. We 
also compute total waiting time of vehicles per cycle behind red lights for various control parameters. 
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Simulation of urban traffic flow has shown to be of 
prime importance for optimisation and control purposes 
as the number of vehicles increases continuously and traf- 
fic conditions deteriorate. Modelling traffic ffow dynam- 
ics by cellular automata has constituted the subject of 
intensive research by statistical physics during the past 
years [ll-[3|- Understanding the characteristics of city 
traffic is one of the most essential parts of the traffic 
research and was an early simulation target for statis- 
tical physicists [341^. Despite elementary attempts for 
simulation of perpendicular ffows Ts'-Tt'] the first serious 
cellular automata describing the ffow at a signalised inter- 
section was proposed by authors in [l8| . Recently, physi- 
cists have notably attempted to simulate traffic flow at 
intersections and other traffic designations such as round- 
abouts [Tol - fs^ . Further recent progresses include various 
features such as self-organised controlling of traffic lights 
and flows in networks [33, 34], the modelling of decision- 
making at intersections 1351 . mixture of motorized and 
non motorized vehicles 13 6| . simulation of traffic flow at 
T-shaped intersections (37l Issj . data collection in city 
traffic [13, application of fuzzy logic in signalisation of 
intersections [40|, visual simulation of vehicular traffic 
[4l| and dynamic route finding strategy [1^. Vehicular 
flow at the intersection of two roads can be controlled via 
two distinctive schemes. In the first scheme, the traffic is 
controlled without traffic lights [H, \^ . In the second, 
signalised traffic lights control the flow. In Ref. [l^, we 
have modelled the traffic fiow at a single intersection with 
open boundary conditions applied to the streets. In a re- 
cent study, a single intersection operating under fixed 
time and traffic responsive schemes has been explored 
with closed boundary condition '3]. In the mentioned 
papers cars were restricted to move straightly and not al- 
lowed to change their directions when they reach to the 
intersection. In this work we incorporate the possibility 
of turning. This incorporation constitutes a crucial step 
for a more realistic description of traffic flow and would 
certainly help us in a proper coordination of traffic lights. 



Consider the traffic flow at the intersection of two 
streets. Each street conducts a uni directional flow and 
has a single lane. Here for simplicity we exclude the pos- 
sibility of overtaking which requires having more than 
one lane and some laws for lane-changing [4^. The flow 
directions are taken south-north in street one and east- 
west in street two. Vehicles can turn when reaching to the 
intersection. A northward moving car can turn left and 
a westward moving car can turn right (towards north) 
when reaching to the intersection. We model each street 
by a chain of L sites. These perpendicular chains inter- 
sect each other at the middle sites ii = «2 = (see Fig.l 
for illustration). The discretisation of space is such that 
each car occupies an integer number of cells denoted by 
Nceii- Cell length is denoted by Ax. The car position 
is denoted by the location of its head cell. Time elapses 
in discrete steps of At and velocities take discrete values 
0, 1, 2, • • • , Vjnax in which v^ax is the maximum veloc- 
ity measured in unit of The probability of turning 
a north- moving car to west is denoted by Psw Corre- 
spondingly the probability of turning a west-moving car 
to north is denoted by pen ■ 

At each step of time, the system is characterized by the 
position and velocity configurations of cars. The system 
evolves under the Nagel-Schreckenberg (NS) dynamics 
[46j . Let us now specify the physical values of our time 
and space units. The length of each car, Lear, is taken 
4.5 metres. Therefore, the spatial grid Ax (cell length) 



equals to 



4.5 



We take the time step At = 1 s. 



Furthermore, we adopt a speed-limit of 74 km/h. In 
addition, each discrete increment of velocity signifies a 
value of Av — ^f'^^ m/s which is also equivalent to the 

T. — 23 



acceleration value. For Nrpii — 5 we have 



cells per time step. Moreover, the acceleration takes the 
comfort value a = 0.9 The value of random braking 
parameter is taken p = O.l throughout this paper. A set 
of traffic lights controls the traffic flow in a flxed-time 
scheme as follows. 
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FIG. 1: Intersection of two uni-directional single-lane streets 
with possibility of car turning at the intersection. 



A. Fixed Time Signalisation of lights 

In this scheme the traffic hghts periodically turn into 
red and green. The period T, hereafter referred to as 
cycle time, is divided into two phases. In the first phase 
with duration Tg, the lights are green for the westward 
street and red for the northward one. In the second phase 
which lasts for T — Tg time steps the lights change their 
colour and become red for the northward and green for 
the westward street. The gap of all cars are updated 
with their leader vehicle except the ones (in each street) 
which are the nearest approaching cars to the intersec- 
tion. These two cars need special attention. For these 
cars gap should be adjusted with the signal in its red 
phase. In this case, the gap is defined as the number 
of cells immediately after the car's head to the intersec- 
tion point -J . If the head position of the approaching car 
lies in the crossing point the gap is zero. Note that for 
Nceii > 1 at a time step when a light goes red for a di- 
rection, portion of a passing car from that direction can 
occupy the crossing point. In this case the leading car 
of the queue in the other direction should wait until the 
passing car from the crossing point completely passes the 
intersection i.e.; its tail cell position become larger than 

L 

2 ■ 



B. Modelling the car turning 

To model car turning we proceed as follows: once a 
car becomes approaching to the intersection we draw 
a random number r (uniformly chosen from the unit 
interval). If r is less than the corresponding turning 
probability the car turning label becomes one otherwise 
it remains zero. In the subsequent time steps the 
approaching car's gap is adjusted appropriately: if its 




(a) (b) 



FIG. 2: A turning adjusts its gap to the appropriate front car 
in the destination lane, (a): the approaching car in street one 
does not intend to turn, (b); the approaching car in street 
one intends to turn. 



turning label is zero its gap is adjusted via its forward 
leading car otherwise if its turning label is one the gap is 
adjusted via the first car in the perpendicular direction 
which has already passed the intersection (see figure 2). 



C. Entrance/Removal of cars to/from the 
intersection 

The type of boundary condition we implement in 
this paper is open, the reason is that if one tries to 
imply close type boundary condition, after some time 
steps most of the cars will appear in that street which 
has a higher turning probability into it. This effect 
arises because of accumulation of cars in the street 
with higher inward turning probability. To avoid such 
spurious effect which is surely contradicting to what 
occurs in reality we impose open boundary condition. 
This type of boundary condition is more compatible to 
realistic traffic flow. For this purpose, we define two 
integer-valued parameters di > 1 and d2 > 1 for streets 
one and two respectively, di is inversely proportional 
to the traffic density in street i. Once the distance of 
the nearest car to first chain cell j — 1 (farthest car to 
the intersection) in street i (i = 1,2) exceeds di a new 
car is inserted at site j = 1 of street i with maximum 
velocity Vmax- The removal of cars from a street is 
modelled as follows: when the position of a car which 
has passed the intersection becomes larger than L — Nceii 
this car is removed from the system and the number 
of exited cars is increased by one correspondingly. The 
output current from street i is defined as the number 
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FIG. 3: Ji versus Tg at Vmax ~ 4 for cycle time T = 
30, 40, 50, 60 and 70 seconds. The values of other parame- 
ters are specified in the figure. Street one has a higher input 
rate than street two. Turning probabilities are equal to each 
other. 

of cars exited from street i during a time interval 
[Ti,T2] divided by the T2 - Ti. We denote the output 
currents by Ji and J2 for street one and two respectively. 



III. SIMULATION RESULTS 

Now all the computational ingredients for simulation 
is at our disposal. We have taken the streets sizes 
equal to 1350 m. For Nceii = 1 this corresponds to 
Li = L2 = 300 cells). The system is updated for 2 x 10^ 
time steps. After transients, two streets maintain steady- 
state currents Ji and J2 which are defined as the number 
of vehicles passing the intersection per time step. They 
are functions of input rates ^ and ^ , turning probabil- 
ities Psw , Pen and signal times T and Tg . 

A. Ncell = 1 

In this case, each cell accommodates a car. Figure 
(3) shows variation of Ji versus green period of street 
two (east-west) Tg for some values of cycle time T. The 
maximum of Ji is less than Tg = ^ because street two is 
less loaded than street one. The inequality of input rates 
makes the slope Ji different at the vicinity of ^. For 
a given set of parameters, one can maximise the output 
current of street one by setting Tg to the maximising 
value. Sharp fall of Ji is noticeable and is related to 
enhancing of jammed region in street one. 

Figure (4) depicts the behaviour of Ji versus Tg but 
this time for various values of input parameter c?i . Other 
parameters are specified in figure. By decreasing the 
green time of street one, its output current slightly in- 
creases until it reaches to plateau region on which Ji is 
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FIG. 4: Ji vs Tg for various values of di at T = 50 and 
Vmax ~ 4. Other parameters values are specified in the figure. 



maintained constant. The reason is that the red cycle 
of light optimally organises the output fiow from the in- 
tersection. This might be seem counterintuitive to the 
expectation the less the red period the larger the out 
flow. However, we should note that turning cars from 
east-west street contribute to Ji too. 

Figure (5) exhibits J2 versus the the green time de- 
voted to street two for various values of input rates of 
street one. Two distinctive behaviour can be identified. 
For a given di the current J2 increases by increasing 
street two green time Tg up to certain value. Afterwards 
J2 remain constant on a small plateau region and start 
to decrease in a slight manner afterwards. The reason 
is that by increasing Tg the green time of street one de- 
creases. This reduces the current of south-west turning 
car into street two. In fact, J2 has two sources: straight 
moving east-west cars in street two and south-west turn- 
ing cars in street one. Beyond a certain Tg the contribu- 
tion of the second source dramatically falls which makes 
J2 show a maximum. 

Let us now look at the behaviour of the total flow 
Jtot — Ji + J2 versus Tg and di. Figure (6) shows this 
variation versus street two green time Tg for various val- 
ues of cycle time T. The notable point is that the maxi- 
mum total flow only weakly depends on the cycle time T. 
For each T there is a plateau region maintaining a max- 
imum flow. This region is centered at a = aT and 
increases with cycle time T. a depends on the streets 
input rates and is less than 0.5 when d2 is less than di. 
Roughly speaking a — j^^- The slope is symmetrical 
on both side of the plateau region. 

In flgure (7) Jtot is sketched versus Tg this time for 
various values of street one input rate di . Typically the 
total flow shows a three-regime behaviour: an increas- 
ing, a plateau and eventually a decreasing portion. The 
plateau length decreases with increase in c?i . These three 
portions are associated to the competition of two compet- 
ing features. The direct current and the turning current. 
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FIG. 5: J2 vs Tg of street two for various values of di at 
Vmax = 4. Other parameters values are specified in the figure. 
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FIG. 6; Jtot vs Tg of street two for various values of T at 
Vmax ~ 4. Other parameters values are specified in the figure. 



The plateau region corresponds to the situation where 
two features are balanced and equally contribute to the 
total current. 

To gain a deeper insight into the problem, we have 
obtained the mutual dependence of the total flow Jtot 
on input rate parameters di and d2 for equal turning 
probability Pgw — Pen = 0.1 in figure (8) (see the jour- 
nal printed version). Other parameters are specified in 
the figure. We see the existence of a two dimensional 
plateau on which the total fiow is maximised. The maxi- 
mum throughput flow is slightly less than 0.8. Moreover, 
the total current Jtot is symmetric respect to the plane 
di = c?2 as it should be. The staircase like topogra- 
phy of the current surface and the appearance of semi 
plateau regions are noticeable. We see the dramatic fall 
of the total current in the vicinity of the plane di — d2. 
This marks the inefficiency of the fixed-time scheme in 
organising the flow via equal distribution of cycle times 
to equally loaded streets. In figures (9) (see the journal 
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FIG. 7: Jtot vs Tg for various values of di at Vmax ~ 4. Other 
parameters values are specified in the figure. 




FIG. 8: Three dimensional plot of Jtot vs T and Tg. di = 
7, d2 = 15, Vmax ~ 4 and turning probabilities are equally set 
at O.f. Jtot retains a constant value on a wide plateau region 
and sharply falls outside the plateau. 



printed version) and (10) we have sketched the Jtot ver- 
sus T and Tg for given input rates. As you can see there 
is a relative large plateau on which the total current is 
maximal. The shape, size and orientation of this plateau 
show strong dependence on input rates c?i and c?2- The 
maximal current value poorly depends on input rates. 

Dependence of Jtot on turning probabilities could serve 
us to achieve a better understanding from the overall flow 
characteristics of our single intersection. Figure (11) ex- 
hibits such behaviour (see the journal printed version). 
The current surface has interesting properties. There is a 
wide plateau region maintaining a large current between 
0.8 and 0.85. Sharp decrease is followed when \psw —Pen\ 
becomes large. To see the reason consider a situation 
where psu, is large whereas Psw is small. Therefore most 
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of the entrant cars to street two move straight forward 
to west. On the other hand most of the entrant cars to 
street one prefer to turn left (towards west) when reach- 
ing to the intersection. This causes a dramatic jam in 
the vicinity of the intersection which leads to a sharp de- 
crease of the total current. In figure (12) (see the journal 
printed version) we have depicted a similar 3D plot this 
times for unequal input rates di ^ d2- Other parameters 
are identical to those in figure (11). The prominent ef- 
fect of having unequal input rates is the reduction of the 
maximal current value in the plateau region from 0.82 
roughly to 0.71. Other characteristics are analogous to 
the case di = d2- Comparison of figure (10) and (11) sug- 
gests that fixed-time signalisation can poorly adjust itself 
to variation of turning probabilities and implementation 
of more efficient adaptive strategies is inevitable. 
Waiting times: 

Despite maximising the output currents is highly de- 
sirable for us but it may not be the ultimate task in 
optimisation of traffic fiow. Another important param- 
eter for optimisation of the traffic flow at intersections 
are aggregate waiting times of cars stopping at queues 
formed in the red cycle of the traffic lights. In our model 
this total waiting time, denoted by WTtot onwards, can 
be simply be evaluated. Once a car reaches to the end 
of a queue the WTtot counter is increased by one at each 
time step. This counting is paused once the correspond- 
ing light goes green. In figure (13) we show the WTtot 
per cycle versus green time of street two for various val- 
ues of cycle time T. The system has evolved for 10^ 
timesteps. The first 20000 timestep has been discarded 
for equilibration. Contrary to total current dependence 
on Tg , here we do not have a minimal plateau region but 
rather there exists a global minimum for WTtot- This 
is helpful to us in achieving the optimisation of traffic 
light by simultaneous maximising the total current and 
minimisation of total waiting time. The results of figure 

(13) and (6) suggest to take the Tgat the end of current 
maximal plateau as the optimal choice for the green time 
distribution. 

To gain a deeper insight, we have drawn a 3D plot of 
the WTtot versus various control parameters in figures 

(14) and (16). In figure (14) (see the journal printed ver- 
sion) the mutual dependence of WTtot on input rates di 
and d^ for T = 30 and Tg — 15 is plotted. The quali- 
tative topography of the WTtot surface resembles much 
with figure (8). For both di and c?2 larger than 13 (light 
traffic in both streets) we have a relatively small WTtot- 
Once either of c?i or c?2 becomes less than 11a dramatic 
increase in WTtot is observed. Figure (15) sheds more 
light onto the problem. This figure is a 2D view of fig- 
ure (14) and clearly shows the staircase-like structure of 
WTtot surface. The sharp fall of WTtot is noticeable. In 
figure (16) (see the journal printed version) total wait- 
ing time per cycle is mutually sketched versus turning 
probabilities for fixed equal input rates. We see there is 
a wide minimal plateau on which the total waiting time 




FIG. 9: Total waiting time per cycle vs Tg. Street one ac- 
commodates a heavier traffic flow (di = 8, (i2 = 10). Vmax ~ 4 
and turning probabilities are equally taken 0.1. 
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FIG. 10: Dependence of WTtot per cycle on di for fixed 
values of d2. Other parameters: Vmax = i, T — 30, Tg = 15 
and turning probabilities are equally taken 0.1. 



retains its constant minimal value. On the boundaries 
where \psw —Pen\ becomes large the waiting time dramat- 
ically begins to increase. The behaviours of total current 
and total waiting times are complementary to each other 
i.e.; the maximum of Jtot coincides with the minimum 
of WTtot- it is our expectation that such complimentary 
between WTtot and Jtot remains valid for the cases in 
which the input rates are not equal. 

We end this section by showing 3D plots of to WTtot 
mutually versus T and Tg. In figure (17) (see the jour- 
nal printed version) WTtot is given for di — d2 — 8- 
There is no minimal plateau sharply separated from high 
values. The curve is very smooth and low value region 
gently crosses over to high values of waiting times. When 
considering non equal input rates, the situation becomes 
worse. In figure (18) (see the journal printed version) 
the WTtot is sketched for = 7 7^ c?2 = 15. Once again 
we see the the region of maximal total current does not 
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FIG. 11: Ji vs Tg at Vmax = 23 for N^m = 5. The values of 
other parameters are identical to figure (3). 
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FIG. 13: J2 vs Tg of street two for various values of di at 
Nceii = 5. Vmax ~ 23. Other parameters values are identical 
to figure (5). 
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FIG. 12: Ji vs Tg for various values of di with N^eii ~ 5. 
Other parameters values are identical to figure (19). 



entirely overlap with the region of minimal waiting time. 
This demonstrates the optimisation of trafhc flow is not 
a simple task at least within the framework of NS model. 
Obviously empirical data is needed and can illuminate 
the problem. 



B. Ncell = 5 

In this section we show the results for a multi-cell oc- 
cupation of cars. Here we take Nceii = 5 as discussed in 
section II. In this case we have the comfort acceleration 
a = 0.9 m/ s^. We have reproduced figures 3-7 to see the 
similarities/differences imposed by varying Nceii ■ Figures 
(19-20) show exactly the same graphs as in figures (3-4). 
The only difference is due to Nceii which is now equal to 
5. 

The most prominent difference is that in Nceii = 5 the 
current Ji is decreasing whereas in N^eii = 1 it slightly 
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FIG. 14: Jtot vs Tg of street two for various values of T at 
Nceii = 5. Other parameters values are identical to figure (6). 



increases with Tg and then starts to decrease. Currents 
Ji in uni-cell i.e., Nceii — 1 are in general larger than 
their counterparts in multi-cell case. This is due to unre- 
aUstic large acceleration in the NS model for uni-cell oc- 
cupation. Figure (21) sketches the same diagram shown 
in figure (5) for Nceii — 5. As you see J2 is an increas- 
ing function of Tg while in uni-call case it reaches to a 
maximum, then after a short plateau it declines. 

Figures (22) and (23) exhibit dependence of total cur- 
rent Jtot versus Tg for various values of cycle time T and 
input rate di. 

Here also we encounter substantial differences. The 
most distinguishable feature is the convex nature of total 
waiting time curve. In sharp contrast to uni-cell case, in 
multi-cell occupation, we do not have a plateau region. 
Instead a smooth decrease is observed up to a minimum 
and then it symmetrically begins to increase. Analo- 
gously, the total current value is less than the uni-cell 
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FIG. 15: Jtot vs Tg for various values of di at Nceii ~ 5. 

due to the same reason explained for Ji and J2. In fig- 
ure (23) and for di = 5 The behaviour for small Tg is 
similar to uni-cell case i.e.; Jtot increases. The difference 
is that in uni-cell it reaches to a maximum and then in- 
crease sharply whereas in multi-cell occupation Jtot does 
not decrease after linear increasing regime is ended. For 
di > 5 we do not have the increasing behaviour for small 

IV. SUMMARY AND CONCLUDING 
REMARKS 

In summary, we have simulated the vehicular traffic 
flow at a single intersection in which the possibility of 
turning when cars reach to the intersection is augmented 
to the problem. A set of traffic lights operating in a fixed 
time scheme controls the flow. We have tried to shed 
some light onto the problem of traffic optimisation by 
extensive simulations. Total current is shown to remain 
constant on a plateau of green time period given to one 



of the directions. Three dimensional plot of total current 
for given input rates of vehicles in terms of signalisation 
parameters shows the existence of a 2D plateau region 
encompassed by almost flat planes of sharp decreasing 
currents. By extensive simulations we have examined 
the effect of turning on the output current. The depen- 
dence of Jtot on the whole range of turning probabilities 
for fixed values of other parameters have been computed 
including the equal and nonequal input rates. The main 
finding is the appearance of a plateau current associated 
the central region of the Pg^ — Pen plane. The current 
sharply decreases when \Psw — Pen \ becomes large. Need- 
less to says the plateau characteristics i.e.; its size and 
height depend on the input rate and signalisation param- 
eters T and Tg. Similar properties is observed if instead 
of Jtot one looks at total waiting time per cycle. Note 
here the waiting times sharply increase when \Psw — Pen\ 
becomes large. Besides total current, total waiting time 
per cycle has been computed. Our investigations reveal 
that in the parameter space, the minimisation of total 
waiting time per cycle does not fully coincide with the 
maximisation of total current. This arises the natural 
question of what quantity should be optimised in order to 
acquire the most efficiency for the intersection ? In our 
Nagel-Schreckenberg cellular automata, not only the case 
of uni-cell occupation, where each cell accommodates a 
car, has been considered but also the more realistic case 
of multi cell occupation in which a car occupies more 
than one cell has been investigated. The correspond- 
ing results differ quantitatively and in some case qualita- 
tively. Seeking for advanced models properly designated 
for modelling the behaviours of drivers at intersections 
is inevitable and unavoidable. We would like to end by 
emphasizing on the role of empirical data for adjusting 
the parameters of any intersection traffic model. Despite 
utilising multi-cell occupation can render the decelera- 
tion/acceleration value of moving cars to a realist one, 
the behaviour of halted cars in the queue when the light 
goes green is still suffering from a satisfactory modelling. 
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